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The latest IceCube data suggest that the all-flavor cosmic neutrino flux may be as large as 
10“^ GeV cm~^ s“^ sr“^ around 30 TeV. We show that, if sources of the TeV-PeV neutrinos are 
transparent to 7 rays with respect to two-photon annihilation, strong tensions with the isotropic 
diffuse 7-ray background measured by Fermi are unavoidable, independently of the production 
mechanism. We further show that, if the IceCube neutrinos have a photohadronic (p7) origin, 
the sources are expected to be opaque to 1-100 GeV 7 rays. With these general multimessenger 
arguments, we find that the latest data suggest a population of cosmic-ray accelerators hidden in 
GeV-TeV 7 rays as a neutrino origin. Searches for x-ray and MeV 7-ray counterparts are encouraged, 
and TeV-PeV neutrinos themselves will serve as special probes of dense source environments. 

PACS numbers: 95.85.Ry, 98.70.Sa, 98.70.Vc 


The astrophTOical high-energy neutrino flux observed 
with IceCube [ll-l2| is consistent with an isotropic dis¬ 
tribution of arrival directions, suggesting a significant 
contribution from extragalactic neutrino sources. Most 
likely, the neutrino signals are generated in the de¬ 
cay of charged pions produced in inelastic hadronuclear 
{pp) and/or photohadronic (py) processes of cosmic rays 
(CRs) Hill. All these processes also predict the genera¬ 
tion of hadronic 7 rays from the production and decay of 
neutral pions. The power of multimessenger constraints 
of astrophysical scenarios has been demonstrated [l^ in 
light of the IceCube and Fermi data [l^. CR reservoirs 
such as starburst galaxies and galaxy clusters or groups 
have been considered as promising sources, and neutrinos 
produced by pp interactions between CRs and gas could 
account for the diffuse flux at > 100 TeV [l^. [l3 - [l^ . 

The contribution of astrophysical neutrinos has been 
studied based on various analysis techniques. By now, 
the strongest significance comes from high-energy start¬ 
ing event (HESE) searches with IceCube QiH- A re¬ 
cent combined likelihood analysis Q sensitive to neutrino 
energies of 10 TeV to 10 PeV suggests the all-flavor flux 
is ~ 10“^ GeV cm“^ s“^ sr“^ around 30 TeV and 

a power-law index Sob = 2.50 ± 0.09 (for oc 7?“®°'’). 
The most recent HESE data also indicate such a soft 
component Q- These observations are consistent with 
an equal contribution of three neutrino flavors [itI - I^ . 

This work considers multimessenger implications of the 
latest IceCube results for an extragalactic origin. As 
shown in Ref. [l^ . the neutrino and 7-ray spectral in¬ 
dex should be s < 2.1-2.2 for a power-law oc 
in contrast to Sob ~ 2.4-2.6. In CR reservoir models ex¬ 
plaining < 100 TeV data, the spectral index should be 
s ^ 2.0 and ^ 100% of the isotropic diffuse 7-ray back¬ 
ground (IGRB) comes from the same neutrino sources. 


challenging the pp scenarios. Our results motivate us to 
study p7 scenarios such as models of choked gamma-ray 
burst (CRB) jets [2l| and active galactic nuclei (AGN) 
cores [22 - 1^ . which are opaque to GeV-TeV 7 rays. 

Connecting v and 7 fluxes.— Hadronic interactions 
of CRs lead to the production of mesons (mostly pions), 
which generates a flux of neutrinos via decay processes 
like 7r+ —?► followed by —>■ e'^VeV^. The neutrino 

energy (in the cosmic reference frame) is related to the 
proton energy Sp as Ev ~ (0.04-0.05)ep. The neutrino 
energy generation rate £uQe„ is given by 
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where SpQsp is the CR generation rate. Here the factor 
3/4 accounts for the 1/4 energy loss for the production of 
in the previous decay chain and K denotes the average 
ratio of charged to neutral pions with K kz 1 for py and 
K EE 2 for pp interactions. The energy-dependent me¬ 
son production efficiency, nhn[l, fpp/p-y], accounts for the 
source environment. The corresponding all-flavor diffuse 
neutrino flux, is calculated as (e.g.. Refs. [^. 1^1 

C f (JiZ 

= ^J (1 -H z)^H{z) 


where E„ is the observed neutrino energy and E[{z) is 
the redshift-dependent Hubble parameter. 

The decay of neutral pions 7r° ^ 27 leads to y-ray 
emission. On production, the neutrino and y-rw energy 
generation rates are conservatively related as (9g 
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where y-ray and neutrino energies are related as ~ 
2ei^. However, the generated 7 rays from the sources 
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FIG. 1: Left panel: All-flavor neutrino (thick blue lines) and isotropic diffuse 7 -ray (thin red lines) fluxes for pp and minimal 
P 7 scenarios of Eqs. @ and © that account for the latest IceCube data from ~ 10 TeV to ~ 2 PeV energies @], where 
s' = Sob = 2.5 is used. While pp scenarios require ej) = 25 TeV with a strong tension with the Fermi IGRB [l3|| . minimal py 
scenarios allow the range of 6-25 TeV (shaded regions) as long as the sources are transparent to 7 rays (see the main text 
for details). Right panel: Same as the left panel, but now showing neutrino fluxes of AGN core and choked jet models from 
Refs. [ 21 I, [^. To illustrate the strength of diffuse 7 -ray constraints, we pretend that the sources were transparent to 7 rays. 


may not be directly observable. First, 7 rays above TeV 
energies initiate electromagnetic cascades in the extra- 
galactic background light (EBL) and cosmic microwave 
background (CMB) as they propagate over cosmic dis¬ 
tances. As a result, high-energy 7 rays are regenerated 
at sub-TeV energies [^. Second, intrasource cascades 
via two-photon annihilation, inverse-Compton scattering, 
and synchrotron radiation processes can prevent direct 7 - 
ray escape [ 2 ^. To see their importance, we temporarily 
assume that the sources are 7 -ray transparent. We will 
see in the following that this hypothesis leads to strong 
tensions with the IGRB, disfavored by the Fermi data. 

In pp scenarios, neutrino and generated 7 -ray spectra 
follow the CR spectrum, assumed to be a power law. In 
CR reservoirs such as galaxies and clusters, a spectral 
break due to CR diffusion is naturally expected [IJ, ll5j . 
Thus, the neutrino spectrum is approximately given by 
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where is the break energy and the softening of the 
spectrum, (5 = s' — s, is expected from the energy depen¬ 
dence of the diffusion tensor [ 9 ^. In pp scenarios, the 
corresponding generated 7 -ray spectrum is also a power 
law into the sub-TeV region [see Eg. ©], where it 
directly contributes to the IGRB jlOfll l and Ref. [lH ob¬ 
tained a limit s < 2 . 1 - 2.2 for generic pp scenarios that 
explain the > 100 TeV neutrino data. The limit is tighter 
(s ~ 2 . 0 ) if one relaxes this condition by shifting e!’, to 
< 30 TeV to account for the lower-energy data . 

Motivated by results of Ref. we calculate the dif¬ 
fuse neutrino spectrum using Eq. o with s = 2 and 
s' = 2.5 and the corresponding 7 -ray spectrum using 


Eq. ([3]). Following Ref. [25|, we numerically solve Boltz¬ 
mann equations to calculate intergalactic cascades, in¬ 
cluding two-photon annihilation, inverse-Compton scat¬ 
tering, and adiabatic losses. In the left panel of Fig. [T] we 
show the resulting all-flavor neutrino and 7 -ray fluxes as 
thick blue and thin red lines, respectively, in comparison 
to the Fermi IGRB and IceCube neutrino data [5j. To 
explain the <100 TeV neutrino data, the contribution to 
the IGRB should be at the level of 100% in the 3 GeV to 
1 TeV range and softer fluxes with s > 2.0 clearly over¬ 
shoot the data. As pointed out by Ref. [l^, this argu¬ 
ment is conservative: the total extragalactic 7 -ray back¬ 
ground is dominated by a subclass of AGN, blazars (e.g., 
Refs. [10, Hll), and their main emission is typically vari¬ 
able and unlikely to be of pp origin [ 10 , HO] . Most of the 
high-energy IGRB is believed to be accounted for by un¬ 
resolved blazars [lO-HOj. Although the IGRB should be 
decomposed with caution, if this blazar interpretation is 
correct, there is little room for CR reservoirs [l0 |. 

In p 7 scenarios, neutrino and 7 -ray spectra depend on 
a target photon spectrum. The effective optical depth 
to photomeson production {fp.y) typically increases with 
CR energy, so that the neutrino spectrum is harder than 
the CR spectrum. However, it cannot be too hard since 
the decay kinematics of pions gives s^Qe^, oc ej as a low- 
energy neutrino spectrum [1^ . In minimal p^ scenarios, 
where neutrinos with s,y < < 25 TeV are produced 

by CRs at the pion production threshold, the neutrino 
spectrum is approximately given by 

SuQe,, OC |^ 2 -s' ~ (minimal py). (5) 

In the left panel of Fig. [2 we show the resulting neu- 
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trino and 7 -ra y sp ectra with the diffuse neutrino flux 
and the IGRB |10l| for a neutrino break in the range 
6-25 TeV. Since the sub-TeV emission is dominated by 
7 rays from cascades in the CMB and EBL, the tension 
with the IGRB can be weaker than in pp scenarios. How¬ 
ever, the IGRB contribution is still at the level of ^ 50% 
for e'l = 25 TeV and reaches ^ 100% for = 6 TeV. 

The spectrum (O can be realized when the target pho¬ 
ton spectrum is a power law with a high-energy cutoff or a 
gray body (see below). We note that specific models have 
larger contributions to the IGRB, by accounting for the 
detailed energy dependence of fpp/p-y, the contribution 
from low-energy GRs, and cooling of charged mesons and 
muons. As examples, we consider hadronic 7 rays in the 
low-luminosity AGN model of Ref. [ 2 ^ (model A), which 
can explain < 100 TeV neutrino data, and the choked 
GRB jet model of Ref. [2l| (model B), although these 
sources are predicted to be opaque to very-high-energy 7 
rays. The right panel of Fig. [1] shows the corresponding 
all-flavor neutrino and generated 7 -ray spectra as thick 
blue and thin red lines. Pretending 7 -ray transparency 
leads to violation of the high-energy IGRB data. 

The limits of the IGRB contribution of pj scenarios are 
expected to become even stronger by identifying addi¬ 
tional point sources or by decomposing the emission into 
contributions from individual source populations. This 
should further constrain the 7 -ray transparent sources 
for e(( = 6-25 TeV, which may still be allowed by the 
Fermi data (cf. left panel of Fig. [T]). On the other hand, 
since the sub-TeV 7 -ray emission comes from cascades 
in the CMB and EBL, the tension with the IGRB can 
easily be relaxed compared to pp scenarios if the sources 
are j-ray dark, i.e. if high-energy 7 rays generated in the 
sources of diffuse neutrinos undergo efficient interactions 
with intrasource radiation. In fact, this is generally the 
case for py scenarios as we will show in the following. 

Connecting py and 77 optical depths.— Let us 
consider a generic source with comoving size r/T (where 
r is the emission radius and L is the bulk Lorentz factor 
of the source). We assume the presence of target photons 
with £( « St/T and spectrum rig/. For rig' oc “ with 
a > 1 , which is valid in most nonthermal objects, meson 
production is dominated by the A-resonance and direct 
pion production. Its efficiency fp^ is given by 

fp~f{Sp) ~ ? 7 p 7 (ck)^p 7 (*"/r)(e(ng')|g'—o.SmpC^e/i/eJj > ( 6 ) 

where ap^ ^ 0.7 x 10“^® cm^ is the attenuation 
cross sectiori (the product of the inelasticity and cross 
section |38l - l^ f. rip^{a) « 2/(1 -|- a), and ea ~ 

0.3 GeV. The energy of protons that typically inter¬ 
act with photons with cosmic reference frame energy 
Et is £p ~ 20 ey ss 0.5r^mpC^eAet~^, leading to St ^ 
20 keV (r/10)^(ej//30 TeV)”^. Thus, the IceCube data 
imply sources with x-ray or MeV j-ray counterparts. If 
target radiation is generated by synchrotron or inverse- 
Gompton emission from thermal or nonthermal electrons. 


low-energy photon spectra can be expressed by power- 
law segments, rigj e'^ “, where a > 2/3 [s^. For rig/^ oc 

e'p®" and a > 1 , the efficiency scales as fp^ oc £p~^, 
and the neutrino spectral index is s = Scr + 1 ~ Q^- For 
(7 < 1 the secondary neutrino and y-ray spectra follow 


the initial GR spectrum, s ~ Scr, above the pion produc¬ 
tion threshold because fp-y becomes approximately con¬ 
stant due to higher resonances and multipion produc¬ 
tion [ 3 ^ [ 4 ^ . A similar scaling is obtained for gray -body 
and monochromatic target photon spectra 33 


Now, in p 7 scenarios, the same target photon field can 
prevent 7 rays from escaping the sources. The relevance 
of two-photon annihilation in GRBs and AGN has been 
considered [il|,[4^. The optical depth to 77 —>■ e+e“ is 


'77 
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where ar — 6.65 x 10“^® cm^ and r]-yy{a) ~ 7/[6a®/®(l -|- 
a)] for 1 < a < 7 [i^, which is the order of 0.1. The 
typical y-ray energy is given by E-y ~ T^m^c^Et~^. 

Neutrino sources considered here include transrela- 
tivistic or relativistic sources like GRBs, pulsars, and 
AGN including blazars. For example, the observed neu¬ 
trino energy is expressed to be = E^/i). -\- z) ~ 
re'y/(l -I- z). Equations ([51) and ([7]) can be used for 
both internal and external photon fields. As shown 
in Refs. [H for reprocessed radiation from ionized 
clouds, the cases of T = 1 are reduced to the formu¬ 
las for external photon fields. Thus, regardless of these 
model details, Eqs. and 0 lead to the following re¬ 
lation I 2 I, m, H : 
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0.01 


( 8 ) 


where e^ is the y-ray energy corresponding to the reso¬ 
nance proton energy satisfying 


2 mgC^ 

- 

rripEA 


GeV 


V 25 TeV/ 


(9) 


Thus, the neutrino data from 25 TeV to 2.8 PeV [^, 
corresponding to the proton energy range from 0.5 to 
^ 60 PeV, can directly constrain the two-photon anni¬ 
hilation optical depth at Ej ^ 1-100 GeV. Importantly, 
Eqs. ([5]) and Q are independent of P and valid for both 
internal and external radiation fields. 

In general, the effective py optical depth fp~y depends 
on source models. But too small values of fpj seem un¬ 
natural since the observed neutrino flux is not far from 
the Waxman-Bahcall 45 JE^ (see also Ref. [ 2 ^) and 
nucleus-survival bounds T7| , corresponding to maximally 
efficient neutrino production in the sources of ultrahigh- 
energy (UHE) GRs. More quantitatively, it is possible 
to obtain general constraints on fpj by comparing the 
observed GR and neutrino fluxes. Recently, Ref. [i^ 
obtained fp~y > 0.01 by requiring that the extragalac- 
tic GR flux does not overshoot the observed all-particle 
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FIG. 2: Neutrino and CR bounds on the optical depth to 
77 —^ e"’"e“ in the sources of diffuse TeV-PeV neutrinos. We 
calculate and fpj as functions of and Sp, respectively, 
imposing fp^ > 0.01. We consider simple power laws with 
a = 2.5 and a = 2/3 for e\, = 6-25 TeV (shaded bands), and 
the gray-body case with the temperature fcT/F^ = 112 eV. 


CR flux « 4 X 10"^ GeV cm-^ at 10 PeV 

(e.g., Ref. [ 43 ). Since the observed CR flux in this en¬ 
ergy range is dominated by heavy nuclei from Galactic 
sources such as supernova remnants, this constraint is 
conservative. The recent KASCADE-Grande data 
suggest that a light CR component may become promi¬ 
nent above the second knee energy at 100 PeV, which 
can be interpreted as the onset of an extragalactic com¬ 
ponent. Using their inferred extragalactic, light CR flux 
Ap'hp Ri 2 X 10“® GeV cm“^ s“^ sr“^ as an upper limit, 
we obtain > 0.1 at £p > 10 PeV |l02| . 

A similar conclusion is drawn by examining nonther- 
mal luminosity densities of known objects. The CR lu¬ 
minosity density of galaxies including starbursts is re¬ 
stricted as SpQep ^ lO^^-lO'^® erg Mpc“^ yr“^ [Ell, 
[ 5 ^ . The luminosity density of x rays {Qx ~ 2 x 
10^® erg Mpc“^ yr“ 1 [13), which are thought to orig¬ 
inate from thermal electrons in hot coronae, can be re¬ 
garded as an upper limit of nonthermal outputs from 
AGN. Adopting SpQsp ^ 2 x 10^® erg Mpc“^ yr“^ as a 
reasonable assumption for CRs from galaxies or AGN, we 
have fp-y > 0.01, independently of the above argument. 

Figure [2] shows comparisons of the effective jry optical 
depth required from the IceCube observation to the cor¬ 
responding optical depth to 77 interactions in the Fermi 
range, related by Eq. (I 8 |). Strictly speaking, Eqs. (H]) and 
(HI) are valid for soft target spectra. To see the robustness 
of our results, following Ref. we perform numerical 
calculations using the detailed cross sections of the two- 
photon annihilation and photomeson production (includ¬ 
ing nonresonant processes). We consider target photon 
spectra leading to = 6-25 TeV (indicated as bands in 
Fig- HI) I which can reproduce minimal pj scenarios. Note 
that adopting lower values of or assuming 7 -ray trans¬ 


parency for models like those shown in the right panel of 
Fig. [1] leads to inconsistency with the Fermi IGRB data. 
The conclusion from Eq. ([ 8 ]) holds even for realistic tar¬ 
get radiation fields, including synchrotron and gray-body 
spectra. 

The high jry efficiency suggested by the IceCube data 
and upper limits on CR luminosity densities suggest that 
the direct 1-100 GeV 7 -ray emission from the sources- 
either leptonic or hadronic-is suppressed. Thus, tensions 
with the IGRB, which are unavoidable for 7 -ray transpar¬ 
ent sources, are largely alleviated or even absent. How¬ 
ever, TeV 7 -ray counterparts could be seen by Cherenkov 
telescopes and the High-Altitude Water Cherenkov Ob¬ 
servatory. For power-law target photon spectra, which 
extend to low energies, is larger than unity beyond 
the Fermi band and as a result the TeV emission from 
the sources should also be suppressed (see Fig. (2). For 
gray-body-like spectra, one could expect point-source 7 - 
ray emission above TeV. The escaping hadronic 7 rays 
are cascaded in the CMB and EBL and could be visi¬ 
ble as extended pair-halo emission in the sub-TeV range 
(e.g.. Refs. (^,123). In this special case, although direct 
point-source emission at 1-100 GeV is still suppressed 
and the tension with the IGRB remains, TeV counter¬ 
part searches can be used as an additional test. 

Summary and implications. — We considered im¬ 
plications of the latest IceCube results in light of the 
multimessenger data. Based on the diffuse v-^ flux con¬ 
nection and CR -7 optical depth connection, we showed 
that the two-photon annihilation optical depth should be 
large as a direct consequence of astrophysical scenarios 
that explain the large flux observed in IceCube. 

There are various implications. Cross correlation of 
neutrinos with Fermi-LAT sources is predicted to be 
weak. Rather, in p 7 scenarios, since target photons are 
expected in the x-ray or MeV 7 -ray range, searches for 
such counterparts are encouraged. Candidate sources of 
hidden CR accelerators include choked GRB jets [Elj and 
supermassive black hole core s [EEl lES ] (see also the 
Supplementary Material [l03l| . which includes Refs, (ssl - 
113), so correlations with energetic supernovae including 
low-power GRBs, flares from supermassive black holes, 
radio-quiet or low-luminosity AGN, and a subclass of 
flat spectrum radio quasars can be used to test the mod¬ 
els. For broadband nonthermal target photon spectra, 7 
rays are suppressed at TeV-PeV as well as 1-100 GeV 
energies. However, if the target photons follow a nar¬ 
row thermal spectrum or are monochromatic in x rays, 
hadronic 7 rays might be seen in the TeV range for nearby 
neutrino sources. Although the obvious multimessenger 
relation between neutrinos and 7 rays no longer exists, 
our findings suggest that cosmic neutrinos play a special 
role in the study of dense source environments that are 
not probed by 7 rays. Larger detectors such as IceCuhe- 
Gen2 [93 sensitive to 10-100 TeV neutrinos would be 
important for the identification of the sources via auto¬ 
correlation of neutrino events [U, [13 . 
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We have assumed that the diffuse neutrino emission 
is isotropic. Even if half of the neutrino flux has a 
Galactic origin, which allows somewhat smaller values of 

~ 2 TeV and fp^, our conclusions remain unchanged. 
Future data on the arrival distribution of starting muon 
events will also be useful. 
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Note added .— After this work came out and was sub¬ 
mitted, Refs. have recently appeared and support 

our argument on pp scenarios. 
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SUPPLEMENTARY MATERIAL 

All the key points are described in the main text, which 
are general and not sensitive to details of the astrophysi- 
cal models. We here describe possible scenarios for TeV- 
PeV neutrino sources that are obscured in GeV-TeV 7 
rays, without going through specific details. For candi¬ 
date sources of CR reservoirs including starburst galax¬ 
ies and galaxy clusters, see Refs, [l^, [Sg and references 
therein. 


Candidates of Hidden Cosmic-Ray Accelerators 

By Equations (1) and (2) in the main text, the diffuse 
(all-flavor) neutrino flux from py sources is estimated to 
be 


~ 0.76 X 10"^ GeV cm"^ s'^ sr" 


X min[l, Jp.y]/sup 




P'^ep 


10"^^ erg Mpc ^ yr 


-3 ,rr-l 


where /sup(< 1 ) is the suppression factor due to meson 
and muon cooling, ^2 is a factor accounting for redshift 
evolution of the source density [d^, . For no redshift 

evolution, we have ^2 — 0.6. For the star-formation his¬ 
tory and flat spectrum radio quasar evolution we obtain 
^2 — 3 and ^2 — 8 , respectively. In the following we will 
discuss specific scenarios in terms of their CR luminosity 
density SpQ^p and photomeson production efficiency fpyy. 

At present, there are several models that can explain 
the 10-100 TeV neutrino data. For a power-law proton 
spectrum, the total CR luminosity density (at z = 0) 
is expressed by Qp = {epQep)'R-p, where TZp{ep) is the 
conversion factor; TZp = ln(e™“/e™™) for Scr = 2 and 

Up = - (e“"Vef‘')^“""]/(scr - 2 ) for 

Scr > 2. In the shock acceleration theory, one typically 
expects ~ FmpC^ or F^WpC^. For example, assuming 
^ gQ pgY^ ^ 2 and = 1 TeV lead to Tip - 
10, while Scr = 2.5 and e™'" = 1 TeV give TZp ~ 10 at 
25 TeV. We hereafter use TZp ^ 10 as a fiducial value, 
although lower e™'" (e.g., ~ 1 GeV) leads to larger TZp. 

Choked jets and newborn pulsars .— Massive star ex¬ 
plosions such as supernovae and GRBs are considered as 
promising sites of CR acceleration. CRB prompt emis¬ 
sion is believed to be high-energy radiation from expand¬ 
ing relativistic outflows launched by a black hole with 
an accretion disk or a fast-rotating magnetar. Parti¬ 
cle acceleration may occur both at internal shocks in¬ 
side a relativistic outflow and a pair of external shocks 
caused by the outflow, via the shock acceleration and/or 
magnetic reconnections [5^. GRBs may explain UHE 
CRs [13, [13, since their integrated 7 -ray luminosity 
density ^ 10 ^"^ erg Mpc“^ yr“^ is comparable to 
the differential UHE CR luminosity density SctQect ^ 
0.5 xlO"^'^ erg Mpc“^ yr~^ at 10^®'® eV. However, stacking 



FIG. 3: All-flavor neutrino fluxes for some 7 -ray obscured 07 
scenarios that may account for the latest IceCube data [a. 
The latest data on upgoing neutrinos are also shown |7l|| . 
We show curves of the AGN core model [13, [zl, choked jet 
model pm, and low-luminosity GRB model [ill.[^. Note 
that model uncertainties are large and contributions to the 
sub-TeV IGRB are sufficiently small in these models. 


analyses for observed GRBs lead to stringent constraints. 
It was shown that classical GRBs can contribute < 1% 
of the observed diffuse neutrino flux [53| . 

However, these limits do not apply to low-luminosity 
GRBs and ultralong GRBs. Low-power GRBs may have 
different origins, and most of them are missed by GRB 
satellites such as Fermi and Swift. Their energy budget 
may be comparable to that of classical GRBs, so it is 
possible that they have a significant contribution to the 
diffuse neutrino flux jp] . Theoretically, lower-power 
jets are more difficult to penetrate the progenitor, so it 
is natural to expect “choked jets” [9^. Although too 
powerful jets lead to radiation-mediated shocks and do 
not allow efficient GR acceleration, since all protons can 
be depleted for meson production, choked GRB jets can 
account for the IceGube data [Il|, Hlj, . Not only 

jets but also newbor n p ulsar winds can serve as hidden 
GR accelerators [hi, [ 6 ^. The pulsar wind with F ~ 10® 
lead to ^ 50 TeV neutrinos in the presence of nonthermal 
target photons generated in the nebula. 

Such jet-driven and pulsar-driven supernovae have 
been suggested as origins of low-luminosity GRBs 
(that are often classified as transrelativistic super¬ 
novae) and hypernovae, whose local rates are ^ 
10^-10^ Gpc"® yr-i jeS] and - 4000 Gpc"® yr'^ [lil, 
respectively. The available energy budget is ^ 4 x 
1046 g]-g Mpc”® yr”^, so we expect SpQ^p ^ 4 x 
10^® 7?.”/ erg Mpc”® yr”^. This does not violate 
the total CR luminosity density of galaxies, SpQcp ~ 
10 ^®- 10 "^® erg Mpc”® yr”^ in the ^ 1-10 GeV range [HI, 
[HU, and it is possible for choked jets and pulsars to 
achieve ^ 10”^ GeV cm”^ s”^ sr”^ (see Fig. [3]). 

In addition, very massive stars born at high redshifts 










lead to black holes and could launch jets (e.g., Ref. 

We here point out that choked jets from such high- 
redshift objects could also give a contribution to the dif¬ 
fuse neutrino flux, as considered in Ref. [toI ]. 

Vicinity of supermassive black holes.- AGN show 
broadband nonthermal emission with major contribu¬ 
tions to the x-ray and 7 -ray backgrounds. About 
10% of AGN is thought to have powerful relativis¬ 
tic jets and these radio-load galaxies have a typical 
density of ~ 10 “^ Mpc“^ and jet luminosity of ^ 
10^^ erg s“^. Their energy budget is estimated to be 


erg Mpc ^ yr ^ 


[73|, implying epQe„ < 


~ 3 X 10^® 

3 X 10^^ erg Mpc“^ yr“^ i- This is comparable to 
the 7 -ray luminosity density of BL Lac objects, Q 7 ^ 
2 X 10"^® erg Mpc“^ yr“^ 24 j. Note that the total CR 


luminosity density is larger by Tip, and large CR load¬ 
ing factors of ~ 10-100 are indeed required to explain 
the observed neutrino and/or UHE CR fluxes [H, It^ . 
The popular jet model explains various multiwavelength 
data of AGN [^, and associated CR acceleration at in¬ 
ner jets, knots, hot spots, and radio bubbles or lobes has 
been discussed for decades (see Ref. [t^ as a recent re¬ 
view). In particular, CR acceleration and neutrino emis¬ 
sions from the inner jets, including effects of both internal 
and external radiation fields, have been studied in detail 
in light of the IceCube data [^, [tI, Izll • However, the 
blazar origin of diffuse neutrinos has already been con¬ 


strained by point-source searches [7|, l71|, |79|, l92| . In addi¬ 
tion, one-zone leptonic and leptohadronic models predict 
very hard spectra and do not explain the 10-100 TeV 
neutrino data [H, [tI, [z 3 ) Izll • 

The situation may be different at deeper regions in 
the vicinity of supermassive black holes. Obviously, the 
photomeson production efficiency is higher at inner radii, 
/p 7 ^ OT ISfl l81|, and x rays can be supplied by the 
black hole accretion disk. The x-ray background orig¬ 
inates from the accretion power of supermassive black 
holes, in particular radio-quiet AGN, and is much higher 
than the 7 -ray background. The 2-10 keV x-ray luminos¬ 
ity density of AGN is Qx ~ 2 x 10^® erg Mpc"^ yr"! [H. 
The X rays are thought to originate mostly from thermal 
electrons in the hot coronae. Although there is no com¬ 
pelling evidence for GR acceleration in such dense envi¬ 
ronments, a fraction of the accretion energy could be used 
for GRs, and radio-quiet AGN and low-luminosity AGN 
can give ~ lO”’^ GeV cm” 


m 2 


given that EpQe^, ^ Qp ^ Qx (see Fig. [3]). Source 
models and CR acceleration mechanisms in the vicin¬ 
ity of black holes are uncertain. First, relativistic jets 
may be common in galaxies hosting supermassive black 
holes, even if they are weaker in radio-quiet and low- 
luminosity AGN (that may include normal galaxies), and 
efficient CR acceleration in compact regions in the jet 
may occur [H, Second, as considered in Ref. [24 |. 
stochastic and/or shear acceleration as well as magnetic 
reconnections may occur in radiatively inefficient accre¬ 
tion flows of low-luminosity AGN disks and radio-quiet 


AGN coronae. Third, in sufficiently low-luminosity ob¬ 
jects starved for plasma, electrostatic acceleration in a 
potential gap formed in the black hole magnetosphere 
may also work [s^. Regarding the luminosity density of 
X rays as an upper limit of nonthermal photon outputs 
from AGN, SpQe^ ^ 2 x 10"^® erg Mpc“^ yr“^ would be 
considered to be a reasonable assumption. 


Multimessenger Limits on Galactic Contributions 

Multimessenger data indicate that the diffuse neutrino 
flux measured in IceCube largely comes from extragalac- 
tic sources. Galactic neutrino emission is excepted to 
be strongly anisotropic except for exotic scenarios like 
emission from the Galactic halo. As shown in Ref. [l^ . 
even these scenarios are constrained by the IGRB: the 
spectral index is required to be s < 2.0 if the Galac¬ 
tic emission is quasiisotropic. In addition, there are up¬ 
per limi ts p laced by CR air-shower arrays such as KAS- 
CADE [SSj j and CASA-MIA [s^ in the TeV-PeV range. 
The isotropic diffuse 7 -ray intensity at E.y ~ 300 TeV 
is limited as < 10“® GeV cm“^ s“^ sr“^. With 

if = 2, an upper limit on the isotropic Galactic halo 
all-flavor neutrino intensity at Ey ^ 150 TeV is esti¬ 
mated to be Ey^y < 1.5 X 10“® GeV cm“^ s“^ sr“^, 
which is < 30% of the all-flavor neutrino intensity at 
Ey = 100 TeV for a spectral index of s = Sob = 2.5. 

A significant contribution may also come from the 
Galactic plane, e.g., by diffuse GRs or unresolved su¬ 
pernova remnants. CASA-MIA [s^ gives Ef^^j < 2 x 
10“® GeV cm“^ s“^ sr“^ at Ej ^ 200 TeV for the re¬ 
gion |5| < 5° and 50° < I < 200°. Assuming the uni¬ 
form neutrino intensity over the entire Galactic plane 
AH, the Galactic plane neutrino intensity is constrained 
as Ey<Py < 2 X 10“® (AH/1 sr) GeV cm“^ s“^ sr“^ at 
Ey ~ 100 TeV, which is only 3% of the all-flavor neu¬ 
trino intensity although some special neutrino sources 
could exist outside the array’s field of view. 

If we consider a neutrino emission region around the 
Galactic center or ridge (e.g., Fermi bubbles), the ob¬ 
servational CASA-MIA limit is weakened to E^^j < 
2 X 10~^ GeV cm“^ s“^ sr“^ at Ej 300 TeV, leading 
to Ey<Py < 3 X 10“® (AH/1 sr) GeV cm“^ s“^ sr“^ at 
Ey ^ 100 TeV. Thus, the Galactic center contribution is 
expected to be < 40-50%. In this case, a stronger upper 
limit (< 25%), which is mostly independent of spectral 
assumptions |88l |. is derived from the distribution of the 
high-energy starting events [§]. 

A separate ht of the IceCube data in the Northern 
and Southern Hemispheres resulted in different best-fit 
power-law indices with Sob 2.0 and Sob 2.56, re¬ 
spectively 0. This could indicate anisotropic emission 
originating from the Milky Way, in particular the Galac¬ 
tic center or ridge. However, this asymmetry is not sig¬ 
nificant at present. Also, muon neutrino limits [s^ and 
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diffuse 7 -ray limits already indicate that Galactic contri¬ 
butions should be less than ~ 25-50%. 

So far, (a) extragalactic CR accelerators such as super¬ 
novae, GRBs, and AGN, (b) extragalactic GR reservoirs 
such as starbursts and galaxy cluster or groups, and (c) 
Galactic sources, have been considered as possible origins 
of the IceCube neutrinos. This work showed that a 7 -ray 
obscured population of CR accelerators is suggested by 
the 10-100 TeV neutrino data. In principle, as hinted 


from the latest IceCube analyses @ , the diffuse neutrino 
flux may come from multiple source populations [TJ, [s^ . 
Although fine tuning is needed, a possible Galactic com¬ 
ponent may alleviate the tension in both pp and p"f sce¬ 
narios. The 7 -ray obscured neutrino source component 
is relevant at < 100 TeV and/or > 100 TeV energies. 
If it is dominant only at low energies, the higher-energy 
hard component may come from CR accelerators or CR 
reservoirs as discussed in Refs. [ill. 11^. 1^ . 


